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, Abstract. Supersymmetric Peccei-Quinn models which provide a suitable candidate for the curva- 

' ton field are studied. These models also solve the /i problem, while generating the Peccei-Quinn 

I scale dynamically. The curvaton is a pseudo Nambu-Goldstone boson corresponding to an angular 

degree of freedom orthogonal to the axion. Its order parameter increases substantially following a 
D ' phase transition during inflation. This results in a drastic amplification of the curvaton perturbations. 

^0 ! Consequently, these models are able to accommodate low-scale inflation with Hubble parameter at 

the Te V scale such as modular inflation. We find that modular inflation with the orthogonal axion as 
curvaton can indeed account for the observations for natural values of the parameters. In particular, 
the spectral index can easily be made adequately lower than unity in accord with the recent data. 
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■ INTRODUCTION 

o : 

^ . The precise cosmological observations of the last decade have established inflation (for 

^ I a review on inflation, see e.g. [1]) as an essential extension of the hot big bang model. 

p However, the case of slow-roll inflation (i.e. the case with inflaton mass <^ H^, the 

p • Hubble parameter at the time when the cosmological scales exit the horizon during 

Qj [ inflation) suffers from the fact that, typically, supergravity (SUGRA) introduces [2, 3] 

' corrections to the inflaton mass of order during inflation. One way to keep the inflaton 

>• ■ mass under control is to use as inflaton a pseudo Nambu-Goldstone boson (PNGB) field, 

^ ! since the flatness of the potential of such a field is protected by a global U( 1 ) symmetry, 

j-j I Such candidates are [4] the string axions, which are the imaginary parts of string moduli 

fields with the flatness of their potential lifted only by (soft) supersymmetry (SUSY) 
breaking. This results in inflaton masses ^ H^. The resulting modular inflation can be of 
the fast-roll type [5]. Fast-roll inflation lasts only a limited number of e-foldings, which, 
however, can be enough to solve the horizon and flatness problems. 

In modular inflation, ~ 1 TeV and, thus, the inflationary energy scale is much lower 
than the grand unified theory (GUT) scale. As a result, the perturbations of the inflaton 
are not sufficiently large to account for the required density perturbations for explaining 
the large scale structure in the universe and the temperature perturbations in the cosmic 
microwave background radiation (CMBR). (For a low-scale inflation model where the 
inflaton perturbations are adequate, see [6].) Thus, a curvaton [7] (see also [8]), i.e. 
another "light" field during inflation, is necessary to provide the observed curvature 
perturbation. However, even the curvaton cannot [9] generically help us to reduce the 
inflationary scale to energies low enough for modular inflation. This is possible only in 
certain curvaton models which amplify [10, 11] additionally the curvaton perturbations. 



We will construct [12] a class of SUSY Peccei-Quinn (PQ) models [13] which possess 
such an amplification mechanism and also solve naturally the strong CP and /i problems. 
We use as curvaton an angular degree of freedom orthogonal to the QCD axion, which 
we will call orthogonal axion (the radial PQ field was used as curvaton in [14]). We 
study the characteristics of the scalar potential in this class of models. We then focus on 
curvaton physics and derive a number of important constraints necessary for the model 
to be a successful curvaton model. Finally, we concentrate on a concrete example of this 
class of models and show that it can indeed work for natural values of parameters. 



IVIODULAR INFLATION 

After gravity mediated soft SUSY breaking, the potential of the inflaton s, which is a 
canonically normalized string axion, is [4] 

V(5) = Vm-^m,V + ---, (1) 

where ^ ('?^3/2^p)^ ^iid nis ~ rn-^/z with m3/2 ~ 1 TeV and Mp ~ 2.44 x 10^^ GeV 
being the gravitino mass and the reduced Planck mass respectively. The ellipsis in (1) 
denotes terms which are expected to stabilize V{s) at * ~ Mp. The inflationary potential 
at the time when the cosmological scales exit the horizon is of intermediate scale: 

yi ~ ^m3/2Mp ~ 10^^-^ GeV (2) 

for which ^ m-^ji- 
In this model, inflation can be of the fast-roll type, where 



s = SiexpiF,AN) with F, = -^^l + --lJ, c=(^^) -1. (3) 

Here, AA^ is the number of the elapsed e-foldings and Si the initial value of the inflaton 
field s. From the above, one can obtain the inflation potential e-foldings before the 
end of inflation as 

y(A^)^y^(i-c-2^^^). (4) 

Even in the fast-roll case, modular inflation keeps the Hubble parameter H rather rigid. 
Indeed, the slow-roll parameter 

because c ~ 1 and s <S Mp during inflation (the overdot denotes derivative with respect 

to the cosmic time 0- 

The initial conditions for the inflaton field are determined by the quantum fluctuations 
which send the field off the top of the potential hill. Hence, we expect that the initial 



value for the inflaton is 



^^^^^ 

Assuming that the final value of s is close to its vacuum expectation value (VEV) 
(5) ~ Mp, we find, from (3), that the total number of e-foldings is 

Mot~^ln(^|, (7) 

where we took into account that //m ^^3/2- 



AMPLIFICATION OF THE CURVATON PERTURBATIONS 

We consider a PNGB curvaton o whose order parameter v = v{t) (determined by the 
values of the radial fields in the model) takes [10, 11] a different (larger) expectation 
value in the vacuum than during inflation and, in particular, when the cosmological 
scales exit the horizon. The potential for the real canonically normalized field a is 



V{o) = {vm^y 



1 — cos 



(7) 



V(|(7| <V) 



1-2 2 



(8) 



where m<7 = ma{v) is the variable mass of a. In the vacuum, v = vq and nia — ma- 
in the curvaton scenario, the curvature perturbation observed by the cosmic mi- 
crowave background explorer (COBE) [15], ^ ~ 2 x 10~^, is given by 



(9) 



where is the partial curvature perturbation of the curvaton and fldec is the ratio of the 
curvaton energy density Pa to the total energy density of the universe p at the time of 
the decay of the curvaton: 



10-2 < Hdec = — 



< 1. 



(10) 



dec 



The lower bound originates (see [16]) from the 95% confidence level bound on the 
possible non-Gaussian component of the curvature perturbation from the CMBR data 
obtained by the Wilkinson microwave anisotropy probe (WMAP) satellite [17]. The 
partial curvature perturbation of the curvaton when the latter oscillates in a quadratic 
potential is given [18] by 



„ 5o 



5a 



dec 



(11) 



where "osc" denotes the onset of curvaton oscillations at a time given by //qsc ~ (we 
assume that m<7 had reached its vacuum value before the onset of oscillations). 



The phase Q = o /v corresponding to a remains frozen until the oscillations begin: 

0osc^0*, 50osc^5e„ (12) 

where star denotes the values of quantities at the time when the cosmological scales exit 
the inflationary horizon and 50 is the perturbation in Q. This implies that 
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(13) 



From the curvaton perturbation during inflation 5(7* = H^/27t, we then find [12] that 

(14) 



5(7osc — -^z — , where e = — < 1. 
27te vq 



So, after the end of inflation when v assumes its vacuum value, the curvaton perturbation 
is amplified by a factor e~^. Finally, one can show [12] that 



Obsc ~ V e > gmin = 



_ ~ —111111 - r\ 1 

TteQ ZTlVQ 

where we used the relations 5 (7* /a* < 1 and (Tosc ~ vq. 



(15) 



CAN WE CONSTRUCT PQ IVIODELS WITH A PNGB CURVATON? 

In the SUGRA extension of the minimal supersymmetric standard model (MSSM), there 
exist D- and F-flat directions in field space which can generate intermediate scales 

Mi ~ (m3/2Mp)^, (16) 

where n is a positive integer. It is natural to identify Mi with the symmetry breaking scale 

/a of the PQ symmetry U(1)pq, such that a fi term is generated with jU ~ /g+^/Mp ~ 
m3/2 [19]. This would simultaneously resolve the strong CP and jU problems of MSSM. 
For this, we need a non-renormalizable superpotential term 

AF'+i/ii/j2/M^, (17) 

where A is a dimensionless parameter, P is a standard model (SM) singlet superfield 
with (P) ~ Ml and h\, h2 are the electroweak Higgs doublets. One can show [12] that P 
must necessarily carry a non-zero PQ charge. As a consequence, P has a flat potential. 
To lift the flatness of its potential and generate an intermediate VEV for P ^ Mi, we 
must introduce [20, 21, 22] a second SM singlet Q with non-zero PQ charge having a 
coupling of the type 

^pn+^-kQk/^n^ (18) 

where is a dimensionless parameter and ^ is a positive integer smaller than n + 3. 



After soft SUSY breaking, the scalar potential possesses [12] non-trivial minima at 

1^1, lel-K/zM^)^, (19) 

where U(1)pq is spontaneously broken and /a and jU are generated dynamically. The 
soft masses-squared mj,, Mq ~ m^^j of Q have either sign, while the coefficient 
A of the soft A-term corresponding to the coupling in (18), which is generally complex 
with I A I ~ m3/2 , rnust be large enough for the non-trivial minima to exist if mj, , Mq > 0. 

To implement our scenario, we need a valley of local minima of the potential which 
has negative inclination and \P\ <^ \P\o and \Q\ <^ \Q\o, where |P|o and \Q\o are the 
vacuum values of l^l and \Q\ respectively. If the system slowly rolls down this valley 
during inflation, the order parameter v <S vq and our amplification mechanism for the 
curvaton perturbations may work. This can be achieved only if one of the masses- 
squared rrip, rriQ is negative. Let us take mj, < and > 0. In this case, the scalar 
potential is [22] unbounded below unless k — I. So, we restrict ourselves to the case 
k=l. One can show [12] that the orthogonal axion in this case acquires a mass of order 
m3/2 during inflation and, thus, does not qualify as a PNGB curvaton. 

The addition of a third SM singlet superfield S, however, with a coupling 

(20) 

where p, q are non-negative integers with p + q<n + 3 and > 3, can drastically change 
[12] the situation allowing the implementation of our mechanism. In the next section, 
we will present a concrete class of models of this category. 



PQ IVIODELS WITH AN AXION-LIKE CURVATON 

We consider a class of extensions of MSSM which are based on the SM gauge group, 
but also possess a global anomalous PQ symmetry U(1)pq, a global non-anomalous 
R symmetry U(1)r, and a discrete Zl" symmetry. Note, in passing, that global con- 
tinuous symmetries can effectively arise [23] from the discrete symmetry groups of 
many compactified string theories (see e.g. [24]). In addition to the usual MSSM super- 
fields hi, h2 (Higgs SU(2)l doublets), (SU(2)l doublet leptons), 4' (SU(2)l singlet 
charged leptons), qt (SU(2)l doublet quarks), and uj, d- (SU(2)l singlet anti-quarks) 
with / = 1, 2, 3 being the family index, the models contain the SM singlet superfields 
P, Q, and S. The charges of the superfields under U(1)pq and U(1)r are 

PQ: ^(-2), Q{2),S{0),h,h2{n + l), 

,n + 3, — 1, ,n + l, 
R : n-^), e(^), hi, h2{0) (21) 

with the "matter" (quark and lepton) superfields having PQ = —{n+l)/2 and R = 
{n+l){n + 3)/4. The integer n is taken to be of the form 



n = 4l+l, where Z = 0, 1, 2,..., 



(22) 



for reasons to be explained below. Finally, under , P changes sign. Baryon (and 
lepton) number is [25] automatically conserved to all orders in perturbation theory as 
a consequence of U(1)r (and U(1)pq). The Z2 subgroup of U(1)pq coincides with the 
matter parity symmetry Z^^, which changes the sign of all matter superfields. 
The most general superpotential compatible with these symmetries is 

W = yeijhhe'j+yuijqMu'j+ydijqihidj 

(n+3)/4 

+XP''+%h2/M^+ £ IkS'^+^-'^^iPQf^/M^, (23) 

k=0 

where yeij, yuij, ydij are the usual Yukawa coupling constants, A, A/t are complex dimen- 
sionless parameters, and summation over the family indices is implied. 



The scalar potential 

The resulting scalar potential for PQ breaking after soft SUSY breaking is 

V = \Fp\^ + \Fq\^ + \Fs\^ + Vsoft, (24) 



where 



and 



are the F-terms, and 



Fp= 2kXk TT-n , (25) 

k=l ^"^p 

(n+3)/4 ^+3-4k^pQ^2k-l p 

Fq= 2kXk — , (26) 

k=i ^^^P 

Fs= in + 3-4k)Xk (27) 

k=0 ^"^P 



Vsoft^ml\P\^ + ml\Q\^ + mj\S\^ + 



(n+3)/4 m+3-4k{pn\2k 

k% ^p 



(28) 



the soft SUSY-breaking terms. Here, the soft SUSY-breaking masses-squared mj,, m^, 
and are of the order of the tny2 and can, in principle, have either sign. However, 
the potential V is [12] bounded below only if mj,, and Mq are positive. For definiteness, 
we will take these two soft masses-squared to be equal, i.e. we will put mj, — Mq = m^. 
Also, for simplicity, we assumed universal soft SUSY-breaking A-terms with |A| ~ m3/2- 
For reasons which will become clear later, we take m| < 0. Therefore, the origin in 
field space (P = 2 = S = 0) is a saddle point of the potential with positive curvature in 
the P and Q directions and negative in the S direction. We will call it trivial saddle point. 




FIGURE 1. Plot of y defined in (24)-(28) in units of M/^^/M^^ with respect to \S\ and \P\, which are in 
units of Mi. We took n = 5, Mi = (mMp)^/^, mj, ~ nig = — m| = m~, A = —9m, and Aq = Ai = A2 = 1. 
Also, \P\ = \Q\ and 0^ = 9p = 9q = (9$, 9p, 9q are the phases of S, P, Q) so that the potential is 
minimized. The trivial and shifted valleys as well as the trivial and non-trivial minima are clearly visible. 



One can show [12] that the potential V has a "trivial" valley of local minima which lies 
on the S axis (i.e. dXP = Q = 0) and is clearly visible in Figure 1 . Moreover, there exists 
a "trivial" minimum on this valley at |5| ~ {m-^i2M^Y/^'^^^\ where U(1)pq is unbroken 
and no jU term is generated. So, we should avoid ending up at this trivial minimum. 
The potential V possesses [12] non-trivial minima too with 

\P\ = \Ql |5|~(m3/2M{5)^, (29) 

where U(1)pq is broken and a /i term is generated. Actually, after soft SUSY breaking 
(included in V), the only global symmetry surviving in a non-trivial minimum is Z^^. 

The shifted valley of minima. We expand V for |5| ^ |P| ~ The leading term is 

^(0) = 1^«+3| +mp\P\ +mQ\Q\ 



«+3 



.2|A||A„„|<^MI_,„s 



(30) 



where Op and Oq are the phases of P and Q respectively and A A(„_|_3) /4 is taken real and 
negative by field rephasing. 

The potential V(o) is minimized with respect to the phases Gp and Gq for 



w+3 



^{Gp + Gq)=0 modulo 2n, \P\ = \Q\, 

Mrr=^- = nS«if |Ar->4(„+2)™^. (31) 



The presence of the term A(„^3-)/4(P2)("+-^)/^/Mp in the superpotential is vital to the 
existence of this "shifted" minimum. In view of Z2 , however, this term can only exist if 
(n + 3)/2 is an even positive integer, which implies the restriction in (22). 

For l^l <S |P| ~ |(2|, the shifted minimum of V(o) is also a minimum of V with respect 
to \P\ and \Q\ at a. practically ^-independent position. Thus, for small values of \S\, we 
obtain a "shifted" valley of minima of V at almost constant values of l^l and \Q\. This 
valley, which is clearly visible in Figure 1, has negative inclination for non-zero and 
small values of due to the negative mass term of S. It starts from the "shifted" saddle 
point which lies at \S\ =0 and |P|, \Q\ equal to their values at the shifted minimum 
of y(Q) . Note, in passing, that a shifted valley of minima was first used in [26] as an 
inflationary path in order to avoid the overproduction of doubly charged [27] magnetic 
monopoles at the end of SUSY hybrid inflation [2, 28] in a Pati-Salam [29] GUT model. 

The PNGB curvaton. The dominant ^-dependent part of V can be expressed as 

= m2|5|2-2|A||A^|ffit^M^cos(^405+^(0P + %)) 



xcos(405-2(0p + 0g)), (32) 

where Qs is the phase of S and AA(„_i)/4 is made real and negative. The potential V^^^ 
is minimized with respect to the phases Op, Oq, Os for 4% + {n — l){6p+ Gq) /2 = 
modulo 2n. Defining the real canonically normalized fields 

(j)p = V2\P\ep, (j)Q = V2\Q\eQ, (j)s = V2\S\ds, (33) 

we find [12] three angular mass eigenstates on the shifted valley for l^l <^ |P| = |(2|: 

1. The axion field a = (<^p — (^g) / V2, which remains massless to all orders. 

2. (ppQ = (0p + ^g)/V2 with mass-squared m^g = \A\{n + 3)^\X^fi_^_^y4\x+/2r^ my2- 

3. 05, our PNGB curvaton (orthogonal axion) o with suppressed mass-squared 



n + —I ^\p\2 



\S 



2 



(34) 

Note that the symmetry is [12] very important for the PNGB nature of 0^. 

SUGRA corrections [2, 3, 30] during inflation can be introduced by simply replacing 
A, m^, and m| by their effective values: 

A=A + caH, n? = n? + cpqH^, fh^s ^ m^g + csH^ = csH^ -\ml\ (35) 



respectively. Here, ca is a complex parameter of order unity, while cpQ and cs are real 
and positive parameters again of order unity. We can arrange the parameters so that 



m| > during the initial stages of inflation. In this case, the shifted saddle point of V 
becomes a local minimum and the system may be initially trapped in it. H decreases 
during inflation and, thus, at some moment of time, this minimum turns into a saddle 
point and the system starts slowly rolling down the shifted valley. During this slow roll- 
over, 05 = <7 is an effectively massless PNGB field which can act as curvaton. 



CURVATON PHYSICS 



The required £. The phase = % ~ 1 corresponding to the curvaton degree of 
freedom remains frozen until the onset of the curvaton oscillations. Hence, we expect to 
have <7osc ~ Ovq. From (15), we then find 



\Mp) ' 



where we have also used that //* ~ m3/2 and vq ~ {m2/2M^yi^'^~^^\ Assuming that o 
decays before big bang nucleosynthesis (BBN), i.e. its decay width T(, ~ 'w^/vq is 
greater than the Hubble parameter //bbn at BBN, one shows [10, 12] that 



i . , 1 , ,5 



where Tbbn ~ 1 MeV is the cosmic temperature at BBN. From (36) and (37), we find 
[12] that 

8 + log(4c/g) 
n > Y , (38) 

7-log(f^L/e) 

which implies that n > 1 for ~ 1. The requirement that F^ > //bbn yields [12] 

n-9 9 + log( ma /TeV) 

m^>10^TeV ^ n< ^ J (39) 

1 - log(m<j/TeV) 

for < 10 TeV. This inequality demands that n < 9 for mcj ^ 1 TeV. So, we get 
1 < n < 9. 

Reheating the universe. We will assume that the curvaton decays after dominating 
the energy density of the universe, i.e. ^dec — 1> which is [12] crucial for avoiding the 
overclosure of the universe by axions (see below). The curvaton dominates [12] when 
ti = Hdom, where ^ 

^dom ~ (^^) min{mc7, Finf} (40) 

with Finf ~ g^f^s /2 being the inflaton decay width (g is the coupling constant of the 
inflaton to its decay products). It can be shown [12] that the requirement that F^ < //dom 



(i.e. o decays after dominating the energy density of the universe) results in the bound 

n-2 

1 /'W3/2\"+^ 3O-logje-21og0 

g>^{^] ^ n>——^^——^ =^ n>2 for 0-1.(41) 
^ e^\Mp J 15 + logg + 21og0 

Hot big bang begins after the decay of the curvaton at a reheat temperature 



7reh~ vrcT^p~m3/2 I 1 >7bbn for n<9. 



(42) 



Diluting the axions. For n > 1, /a f« vq ~ (m3/2M^) > lO^^ GeV. This nor- 

mally leads to axion overproduction overdosing the universe. However, if the curvaton 
dominates the universe before decaying, the entropy generated [31] during its decay 

Rafter . ( H^oA ' ^ ^^'^0 . ^qI ( mA ^ ^^3^ 



r\j I I — rsj 



5'before \ J m^/lMp ^^3/2 

can adequately dilute the axions (see [32]). 

The evolution of v. The order parameter v oc \S\ must be slowly rolling during 
inflation to preserve the approximate scale invariance of the perturbations (see below). 
So, it should follow the equation 

: ■: 0: : V \S\ 1 f I ffl^ I \ 

3HIS\+„4\S\^0 ^ - = \^^=.cs[LA^-i)h. (44) 

Using (4) and the fact that = csH^ ^ ^sH^i^ — e"^^*^"), this equation becomes 

3 Jin 151 e-2^'^x_^-2F,iv 



cs dN i^e-2F,Ar 



(45) 



where //x and A^x correspond to the phase transition which changes the sign of m| during 
inflation. The solution of (45) is 

where = |S|(A^h<) ~ (e/eniin)-f^* and |S|x = |5'|(A^x) ~ H^jlTi from quantum fluctua- 
tions at the phase transition. The contribution to the spectral index tig from the evolution 
ofv during inflation is [12] -//-^(v/v)* < 0. The WMAP bound [17] on Us then implies 

? / < 0.04. (47) 

It is important to note that, in the present case, the negative contribution to ns from the 
variation of v during inflation leads naturally to spectral indices which can be adequately 
smaller than unity in accordance with the recent WMAP results [17]. 



A CONCRETE EXAMPLE 



From (39), and (41) and in view of (22), we see that not many choices for n are allowed. 
In fact, we can only accept the models with n = 5, 9 (i.e. / = 1, 2) with the latter case 
being marginal. Hence, to illustrate the above, we take an example with n = 5 (i.e. / = 1) 
and the curvaton assuming a random value after the phase transition, i.e. ~ 1. 

The bound in (39) suggests that this case is acceptable provided that ma ~ 220 GeV. 
Using (36), we find that £ ~ 10^^'^ (recall that f^dec — 1) and Cmin ~ 

10-12 ^ which 

yields I^Ih, ~ IC^//*- We also estimate [12] A^^ to be about 38 for ^^3/2- The reheat 
temperature turns out to be Treh — 10 MeV, while the entropy production at curvaton 
decay is given by Safter/'^before ~ 10^'^^. From (41), we then conclude that g > 10""^"^. 

The dilution of axions by the entropy produced when the curvaton decays after 
dominating the universe may lead [20] to a cosmological disaster. A sizable fraction 
of the curvaton's decay products consists of sparticles, which eventually turn into stable 
lightest sparticles (LSPs) in models (such as ours) with an unbroken matter parity. The 
freeze-out temperature of the LSPs is much higher than Treh- Thus, the LSPs freeze out 
right after their production and can, subsequently, overdose the universe. This problem 
can be solved [12] by suppressing the Higgsino components of the lighter neutralinos 
and charginos below 1% and taking the curvaton adequately light. 

One can show [12] that all the requirements mentioned above can be satisfied for 

C5,c< 0(10^4) ^ \ms\, m,<0{lO-^)H,. (48) 

The smallness of |m5| is due to the requirement that |5| is slowly rolling during the rele- 
vant part of inflation so that the approximate scale invariance of the density perturbations 
is preserved, whereas the smallness of ntg to the required value of e (» emin)> which de- 
mands substantial variation of |5| from the phase transition during inflation until the 
time when the cosmological scales exit the inflationary horizon. Such a variation can 
be achieved with a large number of e-foldings (A^x ~ 0(10"^) and //tot ~ 0(10^ - 10^)), 
which means that, in our case, modular inflation is not of the fast-roll type. 

CONCLUSIONS 

We constructed SUSY PQ models generating the PQ scale and the /i term dynamically. 
They contain a successful PNGB curvaton whose perturbations are suitably amplified to 
account for the observed curvature perturbations even in low-scale inflationary models 
such as modular inflation where the inflaton is unable to generate these perturbations. 
The spectral index of density perturbations can easily satisfy the recent WMAP bound 
in contrast to other inflationary models. However, due to the very low value of the reheat 
temperature, baryogenesis may be achieved only via some exotic mechanism (see [12]). 
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